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SUMMARY 


This  report  provides  a  prediction  of  the  far  field  overpressure 
encountered  near  the  M-198  155rmn  howitzer  when  fired  with  charge  M203.  The 
purpose  of  this  effort  was  to  determine  the  feasibility  of  simulating  the  far 
field  muzzle  blast  in  terms  of  gas  dynamics  and  to  correlate  the  pressures 
with  field  trial  data.  The  results  are  extremely  encouraging.  The  use  of  gas 
dynamics  to  understand  the  physical  phenaomena  and  subsequent  predictions  of 
pressure  and  energy  fields  In  the  far  field  is  a  natural  extension  applicable 
to  many  weapon  systems. 
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1.  INTRODUCTION 


1.1  BACKGROUND 

As  part  of  the  gun  blast  overpressure  program,  WRAIR  had  collected  and 
analyzed  statistically  a  considerable  number  of  time  pressure  traces  of  the 
blast  field  In  the  neighborhood  of  the  M198.  Despite  the  considerable  effort 
that  had  been  expended,  there  were  still  several  questions  that  hindered  the 
further  progress  of  quantifying  crew  exposure.  Were  the  several  peaks  seen  In 
the  data  due  to  ground  reflections,  multiple  bursts,  reflections  of  the  gun 
body,  or  characteristics  of  the  Instrumentation?  How  could  shot-to-shot  re¬ 
peatability  be  quantified?  How  could  the  Lovelace  shock  tube  experiments  be 
compared  with  the  actual  gun  firings?  Should  dynamic  as  well  as  static  pres¬ 
sures  be  compared?  The  answer  to  these  questions  seemed  to  lie  In  a  better, 
understanding  of  the  physical  phenomena  of  the  blast  field. 

1.2  OBJECTIVE  OF  JAYCOR'S  ANALYTICAL  EFFORT 

Because  JAYCOR  was  already  actively  Involved  In  the  gun  blast  overpres¬ 
sure  program,  both  In  assisting  In  data  collection  and  analysis,  and  because 
JAYCOR  has  a  strong  group  In  fluid  dynamics,  WRAIR  decided  to  pursue  a  limited 
analysis  of  the  phenomena.  Based  on  those  considerations,  JAYCOR  determined 
that  much  of  the  pressure  trace  detail  could  be  understood  In  terms  of  gas 
dynamics  and  the  gun  and  ground  geometry.  The  objective  of  this  Initial  phase 
was  to  determine  the  feasibility  of  simulating  the  far  field  muzzle  blast  and 
to  Interpret  the  field  data  already  taken. 

I 

1.3  TECHNICAL  APPROACH  i 

It  was  decided  to  treat  the  blast  as  an  Independent  spherical  distur¬ 
bance  for  each  propagation  direction  away  from  the  muzzle.  Each  disturbance 
would  be  characterized  by  a  different  energy  deposition,  the  angular  nonuni¬ 
formity  being  due  to  the  action  of  the  muzzle| brake  geometry*  For  the  purposes 
of  demonstrating  feasibility,  the  angular  distribution  was  chosen  by  comparing 
with  experimental  data.  Later,  the  distribution  could  be  found  dlrtrtly  from 
the  charge  characteristics,  ballistics,  and  muzzle  brake  geometry. 
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The  method  of  characteristics  was  chosen  to  integrate  the  time- 
dependent,  gas  dynamics  equations  because  of  its  ability  to  accurately  repre¬ 
sent  the  sharp  wave  front.  Ground  reflections  were  taken  to  be  given  by  per¬ 
fect  superposition,  but  work  has  also  begun  on  an  exact  treatment  of  the  non¬ 
linear  effects. 

The  Initial  effort  was  focussed  on  the  physical  character  of  the  blast, 
but  when  the  results  were  shown  to  be  so  encouraging  the  pressure  prediction 
model  was  coupled  to  a  model  of  the  lung  and  chest  cavity. 

In  addition,  the  various  models,  including  simple  ones  for  the  burn 
chemistry  and  shell  dynamics,  were  combined  into  an  overall  system  description 
with  Input  quantities  that  are  actually  controllable,  such  as  gun  elevation, 
charge  strength,  muzzle  geometry,  and  output  quantities  of  interest  to  the 
overpressure  program,  such  as  pressure  time  histories  and  lung  dynamics.  Many 
of  the  models  are  in  the  earliest  stages  o'  development  but  putting  them  all 
together  allows  the  sensitivity  of  each  to  be  determined  and  gives  a  basic 
tool  into  which  future  Improvements  can  be  added. 

1.4  MAJOR  RESULTS 

The  principal  features  of  the  pressure  traces  can  be  understood  in 
terms  of  Intermediate  strength  gas  dynamics  disturbances,  ground  reflections, 
and  the  geometric  relationship  between  the  gun  and  the  measuring  station.  The 
jaggedness  of  the  measured  signal  Is  not  contained  in  the  present  model  and 
the  duration  of  the  pressure  undershoot  is  not  as  great  as  that  measured.  For 
the  cases  tested,  the  simulation  indicates  the  correct  trend  of  pressure  peaks 
and  timings  when  distance,  angle,  height  above  ground,  and  gun  elevation  are 
varied.  The  qualitative  shapes  of  the  pressure  signals  are  reproduced  through¬ 
out  the  measurement  field. 

The  lung  model  gives  a  way  of  comparing  various  pressure  traces  in 
terms  of  the  internal  dynamics.  The  agreement  seen  between  measured  and  pre¬ 
dicted  pressure  traces  is  repealed  in  the  lung  response.  Thus,  the  differences 
In  pressure  traces  does  not  significantly  change  the  lung  dynamics  adding  con¬ 
fidence  to  the  simulation  method. 

In  short,  several  independent  tests  indicate  that  the  major  physical 
phenomena  are  captured  in  this  approximate  approach. 
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1.5  RECOMMENDATIONS 


The  present  phase  of  the  work  was  only  Intended  to  demonstrate  the  fea¬ 
sibility  of  simulation  In  assisting  the  gun  blast  overpressure  program  and  to 
obtain  preliminary  verification.  It  would  be  natural  to  extend  the  work  by 
subjecting  the  models  to  detailed  verification  and  then  useful  application. 
The  Items  below  are  some  of  the  areas  that  could  profitably  be  pursued.  The 
first  group  concerns  the  simulation  of  the  blaut  field. 

(1)  The  present  effort  has  compared  simulation  with  field  test  data 
for  only  about  20%  of  the  shots  measured.  Procedures  have  been  developed  for 
efficiently  handling  all  of  the  data,  so  that  the  validation  can  be  made  com¬ 
plete. 

(2)  Although  the  Ideal  ground  reflection  approximation  appears  to  be 
adequate,  the  analysis  of  the  nonlinear  effects  already  begun  should  be  fin¬ 
ished  because  the  crew  stands  In  an  area  where  this  Is  most  likely  to  be  1m- 
po**tant. 

(3)  In  the  present  treatment  the  angular  distribution  of  energy  caused 
by  the  muzzle  brake  Is  determined  empirically  from  the  field  data  itself.  BRL 
Is  working  on  detailed  calculations  of  the  flow  within  the  brake  that  can 
eventually  supply  a  more  fundamentally  determined  distribution.  There  will 
always  be  a  need  for  a  simple,  but  quantitative  relation  based  on  the  brake 
geometry,  however,  and  JAYCOR  has  several  approaches  that  could  yield  engi¬ 
neering  quality  results  In  a  relatively  short  time. 

(4)  Investigation  into  the  origin  of  the  apparently  random  pressure 
spikes  should  be  made  to  complete  the  understanding  of  the  physical  phenomena. 
Effects  that  could  be  simulated  and  compared  with  data  Include  randomness  in 
the  source  and  response  of  the  Instrument  stands.  This  noise  component  does 
not  seem  to  affect  the  qualitative  features  of  the  signal  other  than  Its  jag¬ 
gedness,  but  the  apparently  random  superposition  can  occasionally  increase  the 
pressure  peaks  considerably  and  It  would  be  useful  to  quantity  the  probability 
of  occurrence. 
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(5)  The  application  of  the  system  code,  that  is  the  collection  of  en¬ 
gineering  models  from  the  explosion  In  the  powder  chamber,  shell  ballistics, 
brake  redistribution,  far  field  propagation,  and  lung  response,  should  be  pur¬ 
sued  further.  It  Is  possible  that  the  remaining  discrepancies  between  the  pre¬ 
diction  and  data,  mainly  the  shape  of  the  rarefaction,  will  be  removed  by 
solving  the  correct  time-dependent,  boundary- value  problem  for  the  far  field 
rather  than  the  present  Initial  value  problem. 

(6)  The  system  description  could  be  applied  to  other  weapon  systems 
with  suspected  overpressure  problems  to  further  validate  the  concepts  In  those 
cases  where  data  exists  and  to  estimate  the  overpressure  field  where  testing 
has  not  been  done. 

The  second  group  concerns  making  use  of  the  validated  simulation  to 
directly  answer  biomedical  questions. 

(7)  An  accurate  and  complete  description  of  the  blast  field  in  the  vi¬ 
cinity  of  the  Ml 98  and  M109  to  animal  placement,  future  test  instrumentation, 
and  test  environment. 

(8)  A  quantitative  determination  if  blasts  from  shock  tubes  are  equiv¬ 
alent  to  those  of  the  gun  from  the  point  of  view  of  lung  response. 

(9)  An  estimate  of  Internal  organ  dynamics  (such  as  chest  wall  accel¬ 
eration)  to  guide  placement  and  specification  of  Instrumentation  In  animal 
tests. 

(10)  An  Interpretation  of  the  dynamics  measured  under  test  conditions 
In  terms  of  a  fluid-mechanical  model  of  the  organ  so  that  the  results  can  be 
extrapolated  to  other  situations  and  to  man. 
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2.  DISCUSSION  OF  RESULTS 


In  order  to  demonstrate  the  feasibility  of  numerical  simulation,  we 
compare  our  calculated  results  with  those  of  the  measured  gun  blast  overpres¬ 
sure  data  for  M198  Howitzer  and  M-203  charge,  collected  on  November  30,  1978 
and  May  15,  1979  at  Aberdeen  Proving  Grounds.  The  comparison  Is  arranged  In  a 
systematical  way  so  that  the  various  effects,  such  as  the  variation  of  the 
distance  away  from  the  muzzle  brake,  the  angular  variation  of  the  measured 
positions,  the  variation  of  the  height  of  the  measured  positions,  and  the 
variation  of  the  gun  elevation  on  the  pressure  traces  can  be  analyzed.  Subse¬ 
quently,  we  compare  the  results  obtained  from  the  simulation  of  lung  response 
on  the  crew  positions,  based  on  the  calculated  pressure  traces  and  the  mea¬ 
sured  ones,  respectively.  The  method  to  calculate  the  pressure  waves  and  the 
lung  simulation  will  be  described  In  the  appendix. 

A  ground  map  showing  the  measured  positions  Is  given  In  Figure  1.  As  we 
shall  mention  In  Appendix  A,  we  simulate  the  origin  of  the  blast  waves  as  a 
sphere  of  high  pressure  and  high  temperature  gas.  In  order  to  account  for  the 
effect  of  unsymnetrlcal  source  distribution  In  the  muzzle  brake,  we  choose 
different  source  conditions,  e.g.,  strength  and  size  of  the  sphere,  for  each 
radial  line.  The  source  conditions  are  adjusted  so  that  both  the  calculated 
and  the  measured  pressure  traces  are  roughly  matched  at  a  si ngle  position  on 
each  line  and  then  the  calculation  Is  checked  against  the  other  pressure 
traces  at  all  other  positions  along  that  line.  It  should  be  emphasized  that 
the  actual  source  distribution  from  the  muzzle  brake  can  be  obtained  by 
solving  the  problem  of  charge  detonation  and  the  resulting  barrel  and  muzzle 
flow.  In  a  future  study,  we  will  calculate  this  part  of  flow  and  integrate  It 
with  the  far  field  calculation  In  order  to  obtain  a  more  complete  and  detailed 
picture. 

The  pressure  traces  at  different  distances  (25",  30",  35",  40",  50", 
and  60")  along  the  B-llne  with  gun  elevation  at  45°  are  shown  In  Figure  2. 
Those  In  the  left  column  are  reproductions  of  the  measured  data  while  the  re¬ 
sults  from  our  calculation  are  shown  In  the  right  column.  The  experimental 
curves  show  a  certain  amount  of  random  oscillation  whose  origin  is  uncertain 
at  this  time.  However,  a  discernible  secona  peak  due  to  wave  reflection  can  be 
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seen  in  all  those  pictures.  Since  the  path  lengths  of  the  Incident  and  re¬ 
flected  waves,  as  shown  In  Figure  42  in  Appendix  A,  are  different  and  their 
speeds  of  propagation  are  a i most  the  same,  these  two  waves  will  arrive  at  a 
certain  point  at  different  times.  Both  our  calculation  and  the  measurement 
display  this  effect.  From  a  geometrical  point  of  view,  the  difference  of  the 
path  lengths  of  these  two  waves  Is  longer  for  the  points  closer  to  the  muzzle 
brake,  as  it  Is  compared  with  those  further  away  from  the  muzzle  brake.  As  a 
consequence,  the  timing  between  the  two  peaks,  as  shown  In  Figures  2  and  3,  Is 
the  largest  at  the  position  with  distance  of  25"  while  It  Is  the  smallest  at 
60".  Due  to  the  longer  path  length  of  the  reflected  wave  than  that  of  the  In¬ 
cident  wave.  Its  strength  should  be  weaker  than  that  of  the  Incident  one.  A 
question  may  arise  as  to  why  the  second  peak  has  a  larger  amplitude  than  the 
first  one  for  some  cases,  for  Instance,  50"  and  60".  This  can  be  explained  as 
due  to  the  overlapping  of  the  resident  times  for  both  waves,  their  amplitudes 
are  superimposed  to  attain  such  a  value,  exceeding  that  of  the  first  one.  Gen¬ 
erally,  our  calculation  has  been  quite  successful  in  capturing  the  structure 
of  the  compression  part  of  the  wave.  As  shown  In  Figure  3,  good  agreement  can 
be  seen  among  the  calculated  and  measured  quantities,  e.g.,  the  amplitude  of 
the  Incident  wave  pi,  the  ratio  of  the  amplitudes  of  the  second  peak  to  the 
first  one  p2/pl,  and  the  timing  between  the  two  peaks  at.  Nevertheless,  the 
rarefaction  part  of  the  wave  In  our  calculation  Is  narrower  in  time  and  shal¬ 
lower  In  amplitude  than  that  of  the  measured  one.  This  discrepancy  Is  probably 
due  to  our  Inexact  treatment  of  the  time  dependent  boundary  condition  at  the 
muzzle  brake,  which  can  only  be  obtained  after  the  abovementl oned  problem  of 
charge  detonation  and  Its  resulting  barrel  flow  Is  solved. 

Figures  4  and  5  show  the  lung  response.  In  terms  of  lung  pressure  and 
Its  change  In  gaseous  volume,  on  the  crew  positions  along  the  B-llne  at  45# 
gun  elevation.  Those  In  the  left  column  are  the  results  from  the  lung  simula¬ 
tion  with  the  measured  overpressure  as  its  driving  force  while  those  using  our 
calculated  overpressure  as  driving  forces  are  shown  In  the  -Ight  column.  Very 
good  agreement  can  be  seen  again  In  the  compression  part.  In  response  to  the 
random  oscillation  associated  with  the  measured  pressure  traces,  as  mentioned 
above,  the  corresponding  lung  pressure  and  Its  change  of  gaseous  volume  show 
wiggles  In  both  the  expansion  and  tall  regions.  On  the  other  hand,  a  regular 


Figure  5.  Change  of  gaseous  volume  of  lung  based  on  measured  and  calculated 
blast  pressures  along  B-llne  at  45*  gun  elevation 
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underdamping  In  the  tall  region  can  be  seen  In  the  counterparts  with  our  cal¬ 
culated  pressure  traces  as  the  driving  forces.  It  should  be  noted  that  the 
validity  of  the  lung  model  Is  established  based  on  the  sound  comparison  be¬ 
tween  numerical  simulation  and  experiments  on  animals  [Reference  1],  Since 
there  Is  no  experiment  on  human  beings  to  our  knowledge,  It  seems  logical  to 
assume  the  validity  of  this  model  can  be  extended  to  study  the  response  on 
human  lung.  In  light  of  this,  let  us  suppose  those  In  the  left  column  are  the 
true  human  lung  response  to  the  real  blast  overpressure,  then  our  calculation, 
as  shown  In  the  right  column.  Indeed  predicts  the  right  trend. 

Figures  6  and  7  show  the  pressure  traces,  pi,  p2/pl  and  At  along  the  C- 
llne  at  45*  gun  elevation.  Figures  8  and  9  show  the  corresponding  lung  pres¬ 
sure  and  change  In  gaseous  volume.  Good  agreement  can  be  seen  again.  The  re¬ 
producibility  of  the  measured  data  and  Its  corresponding  lung  response  on  crew 
position  at  C22  with  45°  gun  elevation  Is  shown  In  Figure  10  for  four  differ¬ 
ent  shots.  The  general  structure  of  these  curves  Is  well  preserved,  except  no 
distinguishable  second  peak  can  be  seen  In  the  pressure  trace  for  Shot  No.  2. 
Moreover,  Figure  11  shows  very  close  agreement  between  the  measured  and  calcu¬ 
lated  pressure  traces  and  their  corresponding  lung  responses  on  crew  position 
at  D60  with  45°  gun  elevation. 

The  same  sequence  of  comparison  from  Figures  2  to  11  for  45°  gun  eleva¬ 
tions,  as  mentioned  above.  Is  repeated  for  15"  and  Is  shown  In  Figures  12  to 
21.  In  view  of  the  gun  geometry,  the  difference  of  the  path  lengths  between 
the  Incident  and  reflected  waves  Is  shorter  for  modest  gun  elevation,  as  It  is 
compared  with  that  of  higher  elevation.  Consequently,  the  timings  between  two 
peaks  are  shorter  for  15°  gun  elevation  than  those  with  45*  (see  Figures  13 
and  17  versus  Figures  '  and  7  for  comparison).  It  Is  the  same  reason,  together 
with  the  argument  of  the  overlapping  of  the  resident  times  of  the  two  waves, 
that  explains  why  the  second  peak  overtakes  the  first  one  for  positions  closer 
to  the  muzzle  brake  at  modest  gun  elevation.  Again,  Figure  20  shows  the  repro¬ 
ducibility  of  the  measured  data  and  Its  corresponding  lung  response  at  C22 
with  15°  gun  elevation  for  four  different  shots.  These  curves  look  very  simi¬ 
lar  with  an  exception  that  the  pressure  trace  In  5hot  No.  23  has  smaller-  sec¬ 
ond  peak  while  the  othe'*  three  cases  do  not. 
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Figure  9.  Change  of  gaseous  volume  of  lung  based  on  measured  and  calculated 
blast  pressures  along  C-llne  at  45°  gun  elevation 
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Figure  14.  Lung  pressure  traces  based  on  measured  and  calculated  blast 
pressures  along  B-line  at  15°  gun  elevation 

21 
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So  far  we  have  shown  In  the  crew  positions  (5-  height)  various  effects 
on  pressure  traces  and  their  corresponding  lung  responses,  namely,  the  varla- 
tlon  of  the  distance  away  from  the  muzzle  brake,  the  angular  variation  of  the 
measured  positions,  and  the  variation  of  the  gun  elevation.  In  order  to  com¬ 
plete  the  study  at  the  crew  positions,  we  Include  the  effect  of  the  variation 
of  the  height  of  the  measured  positions.  Figures  22  and  23  show  the  pressure 
traces,  pi,  p2/pl  and  at  at  different  heights  (3",  4",  5",  and  6“)  at  C22  with 
15*  gun  elevation.  Good  correlation  among  measurement  and  calculation  can 
again  be  observed.  Since  the  difference  of  the  path  lengths  between  the  Inci¬ 
dent  and  reflected  waves  Is  shorter  for  the  points  closer  to  the  ground,  the 
timing  between  the  two  peaks,  as  shown  In  Figure  23,  Is  the  smallest  at  the 
position  with  3"  height  while  It  Is  the  largest  at  6'. 

Figures  24  to  35  show  the  pressure  traces,  pi,  p2/pl  and  at  at  differ¬ 
ent  distances  (10  m,  20  m,  30  m  and  40  m)  along  the  0°,  30°,  60°,  90°,  120® 
and  150®  lines  at  45®  gun  elevation.  In  these  figures,  a  distinguishing  fea¬ 
ture  is  shown,  namely,  the  rarefaction  part  of  the  wave  Is  flattened  out  as  It 
moves  further  away  from  the  muzzle  brake.  The  reason  Is  probably  that  at  large 
distances  the  wave  front  has  a  smaller  curvature  and  behaves  more  like  a  plane 
wave.  Theoretically,  the  plane  wave  consists  only  of  either  a  compression 
pulse  or  a  rarefaction  pulse.  Therefore,  the  waveform  at  points  very  far  from 
the  muzzle  brake,  resembles  the  plane  wave  rather  than  the  spherical  wave. 
Furthermore,  the  timing  between  the  two  peaks  of  the  pressure  traces  becomes 
less  distinguishable  for  points  further  away  from  the  muzzle  brake  since  the 
path  lengths  of  both  the  Incident  and  reflected  waves  are  almost  Identical  at 
those  points. 

Finally,  Figure  36  shows  the  comparison  of  the  calculated  and  measured 
amplitude  of  the  Incident  wave,  l.e.,  the  first  peak  In  the  pressure  trace.  In 
\  the  whole  field  at  45®  gun  elevation.  Both  results  Indicate  that  the  source 

distribution  In  the  muzzle  brake  Is  stronger  near  0®  and  also  In  the  range 
from  60®  to  120®,  which  Is  consistent  with  the  geometrical  configuration  of 
the  muzzle  brake. 
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Figure  25.  Variation  of  pi,  p2/pl  and  At  with  distance  along  0°-line  with  45 
gun  elevation 
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Figure  28.  Measured  and  calculated  pressure  traces  along  60#-11ne  at  45°  gun 
elevation 
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Figure  33.  Variation  of  pi,  p2/pl  and  at  with  distance  along  12G#-11ne 
45°  gun  elevation 
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APPENDIX  A 


A.  GUN  BLAST  FAR  FIELD  ANALYSIS 


A.l  INTRODUCTION 

A  gun  blast  wave  is  a  physical  phenomenon  that  occurs  when  the  high 
pressure  and  high  temperature  product  gas  leaves  the  barrel  after  the  explo¬ 
sion  of  the  charge.  The  atmosphere  is  disturbed  during  and  after  the  emergence 
of  this  gas.  For  the  sake  of  clarity,  we  first  consider  the  blast  waves  in  a 
uniform  unlimited  medium,  vrfiich  can  be  treated  as  spherical  since  they  result 
from  release  of  a  large  amount  of  energy  from  a  source  of  relatively  small  di¬ 
mensions.  The  effects  due  to  the  wave  reflection  from  the  ground,  the  geome¬ 
tries  of  the  gun  and  the  muzzle  brake,  and  so  forth,  will  be  considered  In  the 
following  section. 

At  the  front  of  the  blast  wave,  the  presssure  p  and  density  p  jump 
abruptly  from  their  undisturbed  values.  Immediately  after  the  front  passes, 
this  disturbance  decays  to  zero  very  fast  and  is  followed  by  a  rarefaction 
wave  (Figure  37).  The  structure  of  the  blast  wave,  e.g.,  the  amplitude  and  the 
duration  of  both  the  compression  and  the  rarefaction  parts,  varies  according 
to  its  source  strength  arid  the  ambient  condition.  Suppose  the  ambient  condi¬ 
tion  is  the  standard  atmosphere,  then  an  extremely  strong  source,  which  re¬ 
sults  from,  say,  a  nuclear  explosion,  will  generate  a  strong  blast  wave.  The 


Figure  37.  Spatial  Pressure  Distribution 
in  81ast  Wave 
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sound  waves  are  essentially  constant.  Mathematically,  it  Is  much  easier  to 
treat  the  sound  waves  since  the  governing  equations  can  be  linearized  by  drop¬ 
ping  higher  order  terms.  Also;  they  reduce  to  a  simple  wave  equation.  The 
technique  for  solving  this  equation  can  be  found  in  some  of  the  standard  text¬ 
books  {for  example,  Ref.  5). 

Although  human  beings  have  used  the  gun  powder  for  centuries,  the  blast 
wave  theory  has  not  been  studied  until  World  War  II.  Since  then,  a  large 
amount  of  experimental  as  well  as  theoretical  work  has  been  done.  Perhaps,  the 
greatest  progress  was  made  when  the  concept  of  similarity  was  introduced  inde¬ 
pendently  by  Sedov  [1946,  Ref.  6]  and  Taylor  [1950,  Ref.  7],  The  concept  has 
been  used  In  other  branches  of  fluid  dynamics,  such  as  boundary  layer  theory, 
conical  flow  theory,  and  transonic  and  hypersonic  flow  theory.  The  assumption 
of  similarity  decreases  the  number  of  Independent  variables  and  thus  often  re¬ 
duces  the  governing  partial  differential  equations  to  more  manageable  ordinary 
differential  equations.  It  should  be  noted  that  similarity  in  the  blast  wave 
phenomenon  does  not  hold  exactly,  but  Is  valid  only  when  the  wave  is  strong 
enough  to  neglect  the  effect  of  the  ambient  atmospheric  pressure.  Neverthe¬ 
less,  the  theoretical  study  of  the  blast  wave  phenomenon  has  to  face  the  very 
difficult  problem  of  finding  the  solution  to  the  unsteady  flow  of  the  fully 
nonlinear  gasdynamic  equations,  satisfying  a  moving  boundary  condition  at  the 
wave  front.  Analytical  solution  can  be  obtained  only  for  very  strong  blast 
waves  for  which  the  similarity  assumption  can  be  justified. 

For  blast  waves  with  Intermediate  source  strength,  numerical  methods 
have  to  be  employed  for  solving  the  gasdynamlcs  equations.  Probably  one  of  the 
best  methods  Is  the  method  of  characteristics  which  is  applicable  to  hyper¬ 
bolic  type  partial  differential  equations.  We  do  not  intend  to  go  into  the 
mathematical  theory  of  hyperbolic  equations,  which  may  be  found  in  several 
excellent  works  (for  example.  Ref.  8).  However,  we  shall  simply  mention  the 
main  results  needed  for  the  computation.  The  distinguishing  property  of  the 
hyperbolic  equations  Is  the  existence  of  certain  characteristic  directions  or 
lines  in  the  r-t  plane,  usually  called  characteristics.  It  must  be  noted  that 
the  characteristic  network  Is  not  known  a  priori,  a  result  of  the  nonlinearity 
of  the  equations.  Along  the  characteristics,  the  dependent  variables  satisfy  a 
certain  relation  known  as  the  compatibility  relation.  It  provides  the  key  to 
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the  method  of  computation.  Detailed  procedures  will  be  given  In  the  following 
section. 

In  light  of  the  experimental  gun  overpressure  data  collected  for  M198 
Howitzer  and  M-203  charge,  for  example,  see  Figure  2,  and  the  pressure  traces 
shown  In  Figure  38,  one  can  Immediately  Identify  that  the  gun  blast  waves  are 
blast  waves  with  Intermediate  source  strength.  It  can  be  shown  that,  for  the 
strong  blast  wave,  the  attenuation  of  the  maximum  pressure  at  the  front 
behaves  like  r"2  while  for  the  sound  wave  It  behaves  like  r“*  [Ref.  53.  A 
systematic  analysis  of  the  experimental  gun  overpressure  data  Indicates  that 
the  maximum  overpressure  attenuates  like  r-1  to  r"2.*  Therefore,  the  validity 
of  the  assumption  of  the  gun  blast  waves  as  blast  waves  with  Intermediate 
source  strength  Is  again  Indicated. 

A.2  MATHEMATICAL  FORMULATION 

After  the  explosion  of  the  charge,  the  shell  Is  propelled  by  the  high 
pressure  and  high  temperature  product  gas,  which  leaves  the  muzzle  brake  as 
the  shell  Is  launched.  Initially,  the  gas  flow  In  a  small  region  surrounding 
the  muzzle  brake  Is  complicated  by  the  geometries  of  the  gun  and  the  muzzle 
brake.  This  effect,  however,  becomes  less  significant  further  away.  Suppose  we 
choose  the  center  of  the  muzzle  brake  as  the  origin,  we  can  construct  a  spher¬ 
ical  coordinates  (r,  a,  *)  as  shown  In  Figure  39. 


*For  Instance,  the  following  table  shows  the  measured  values  of  the  amplitude 
of  the  Incident  wave  obtained  from  the  left  column  In.  Figure  32  and  Its  com¬ 
parison  with  the  values  obtained  from  the  ~r“*  and  ~r“2  relationships. _ 

Distance  Distance  from  1  1 


Oeflned  In 
Figure  1 
(m) 

the  Muzzle 
Brake,  r 
(m) 

~r-l 
—  (ps*g) 

1 

1 

Experimental  ~r-2 

_ (pslg)  ~ _ (pslg) 

10 

11.44 

1.41 

2.50  1 

4.98 

20 

20.76 

0.78 

1.06 

1.51 

30 

30.51 

0.53 

0.71 

0.71 

40 

40.38 

0.40 

0.40 

0.40 
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Figure  39.  Spherical  Coordinates  with  Origin  at  Center 
of  Muzzle  Brake 


Assuming  the  gas  Is  Ideal,  Invlscld  and  It  flows  Isentroplcally,*  the 
equations  of  continuity,  motion  and  energy  can  be  written: 

at  +  ^ar1"^  *  FsTrTe  la  {fiV  sin  9)  +  r”sTn“e  I?  M  *  0  (1) 

fbu  v2  +  w2!  _  3P  /-% 


where 


IOu  v2  +  w2]  _  3p 
Bt - r  J  ‘  3? 

SDv  .  uv  w2  cot  el  _  1  ap 

Dt  +  r - r  J  *  "  “r  le 

IDw  .  uw  .  vw  cot  el  _  1  3p 

w  +  7_  +  _ - j  -  *  r”"sTn"""e  "5^ 


Bt  Cp3  ■  *2  W  M 


0  -i_  +  w  i_ 

T5t  a  at  ar  r  ae  r  sin  e  av  * 


♦Isentroplc  assumption  Is  not  valid  for  the  shock  wave  associated  with  the 
strong  blast  for  which  the  Ranklne-Hugonlot  shock  relations  have  to  be 
employed  to  account  for  the  entropy  jump,  but  It  Is  valid  here. 
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Here  u,  v  and  w  are  the  velocity  components,  p  Is  the  pressure,  p  Is  the  den¬ 
sity,  a  Is  the  local  sonic  speed,  and  t  Is  the  time  measured  from  the  Instant 
that  the  blast  waves  are  generated. 

A  complete  solution  to  Equations  (1)  to  (5)  depends  on  the  source  dis¬ 
tribution  In  the  muzzle  brake  which  In  turn  depends  on  the  detonation  of  the 
explosive  charge  and  Its  resulting  flow  field  In  the  barrel.  However,  for  the 
present  study,  our  primary  goal  Is  to  establish  the  feasibility  of  obtaining 
the  solution  of  the  pressure  traces  at  different  far  field  positions  by 
solving  the  gasdynamlcs  equations.  Instead  of  attempting  to  solve  the  complex 
barrel  flow,  we  simulate  the  origin  of  the  blast  waves  as  a  sphere  of  pressur¬ 
ized  gas.  Furthermore,  we  assume  that  the  transport  processes  are  Important 
only  along  the  radial  direction  and  those  associated  with  e  and  41  directions 
are  secondary.  This  assumption  Is  justified  by  the  sound  comparison  of  our 
results  with  the  experimental  data,  as  given  In  Section  2.  With  these 
assumptions.  Equations  (1)  to  (5)  can  be  simplified  as  follows: 

pt  +  Upr  +  p(ur  +  ^r)  “  0  »  (6) 

ut  +  uur  +  J  Pr  *  0  »  (7) 

pt  +  upr  -  a2(pt  +  upr)  «  0  ,  (8) 

where  the  subscript  t  and  r  denote  partial  differentiation.  Here  Equations  (3) 
and  (4)  are  neglected  since  these  equations  describe  secondary  phenomena. 

The  Initial  and  boundary  conditions  to  Equations  (6)  and  (8)  can  be 
written  as 


for  r  >  r0 


(10) 


where  rQ  ■  Initial  radius  of  the  pressurized  gas  sphere,  and  the  subscripts  s 
and  o  denote  source  and  undisturbed  conditions,  respectively. 

At  t  >  0+, 

u  •  0  at  r  ■  0  (11) 

Since  Equations  (6)  to  (8)  are  hyperbolic  equations  [Ref.  8],  their 
characteristic  form  can  be  obtained  In  the  following  way.  First,  let  us  define 

0  3  J  ” f  •  (12) 

Substitution  of  Equation  (12)  Into  Equations  (6)  to  (7)  yields 

ot  +  uor  +  aur  +  -  0  (13) 

ut  +  uur  +  aor  *  0  •  (14) 


Adding  and  subtracting  Equation  (13)  with  Equation  (14),  we  obtain 

(o  ♦  u)t  +  (u  +  a)  (c  +  u)r  +  ~  ■  0  , 

(cr  -  u)t  +  (u  -  a)  (o  -  u)r  +  ■  0  , 

respectively.  It  may  be  shown  that  for  Isentroplc  flow 
2a 


Equations  (15),  (16)  and  (8),  with  the  substitution  of  Eq. 
then  be  cast  Into  characteristic  form  as  follows 
On  r+  curve: 


dr 

US  ”u  +  a 


d  _ 

(IF  T 


— — +  u  ^ 
-  1  r 


0 


(15) 

(16) 


(17) 
(17),  can 


(18) 

(19) 


On  r"  curve: 


I  dr 
3t 

d 

Fr 


■  u  -  a 


2a 


r* 


u 


0 


(20) 

(21) 


On  r°  curve: 

($•« 


(22) 

(23) 


Equations  (18),  (20)  and  (22)  define  the  direction  of  the  character¬ 
istics  r*.  r"  and  r°  respectively.  Along  these  characteristics,  the  compati¬ 
bility  relations,  e.g..  Equations  (19),  (21)  and  (23),  are  satisfied  accord¬ 
ingly.  A  sketch  of  r+,  r"  and  r°  characteristics  are  shown  In  Figure  40.  It 
should  be  noted  that  the  characteristics  r®  are  Identical  to  the  streak  lines 
of  the  fluid  particle.  Furthermore,  the  blast  wave  front  follows  closely  with 
one  of  the  r+  characteristics. 
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Figure  40.  Characteristics  In  r-t  Plane. 


Since  Equations  (18)  to  (21)  consist  of  only  two  unknown  variables  u 
and  a,  these  equations  are  solved  first.  Once  u  and  a  are  found.  Equations 
(22)  and  (23)  can  be  used  to  obtain  the  solution  of  p  after  the  substitution 
of  p  *  vp/a2*  The  solution  of  the  characteristic  equations,  Eqs.  (18)  to  (21), 
Is  based  on  the  upstream  Interpolation  scheme  Introduced  by  Belotserkovskil 
and  Chushkln  [Ref.  9].  The  basic  Idea  Is  as  follows:  Firstly,  approximate 
those  equations  by  use  of  first  order  Implicit  finite  difference  formula  for 
small  time  steps,  which  results  In  a  piecewise  linear  characteristic  network. 
The  wave  front  Is  traced  by  following  one  of  the  r+  characteristics.  Initiated 
at  the  boundary  of  the  pressurized  gas  sphere.  This  particular  characterlst.c 
divides  the  disturbed  flow  region  from  the  undisturbed  one.  A  number  of  points 
In  r-dlrectlon  with  equal  spacing  ar  are  used.  The  exact  number  depends  on  how 
far  from  the  muzzle  brake  we  want  to  calculate.  Old  convergent  solutions  are 
extrapolated  to  give  Initial  guesses  for  flow  variables  at  new  time  level. 
They  are  substituted  Into  Equations  (18)  and  (20)  which  give  roughly  the  loca¬ 
tions  of  the  r+  and  r*  curves  at  old  time  level.  Hence  flow  variables  at  these 
locations  can  be  determine^  by  Interpolation  between  the  convergent  solutions 
at  the  old  time  level.  With  this  Information  fresh  values  at  new  time  loca¬ 
tions  are  obtained  from  Equations  (19)  and  (21).  Due  to  the  Implicit  nature  of 
the  method,  the  process  Is  iterated  until  two  successive  guesses  of  the  flow 
variables  at  the  new  time  Step  agree  within  sufficiently  close  limits.  The 
solution  procedure  can  be  summarized  In  the  following  flow  chart  (Figure  41). 
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**Tne  superscript  v  denotes  the  vth  iteration,  and  c  Is  the  tolerance  limit  of 
relative  error  for  convergence. 

Figure  41.  Flow  Chart  for  the  Method  of  Characteristics. 
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The  wave  reflection  from  the  ground  will  generate  another  peak  In  the 
pressure  traces  (see  Figure  2)  due  to  different  times  of  arrival  of  the  Inci¬ 
dent  wave  and  the  reflected  wave.  The  way  to  calculate  the  reflected  wave  Is 
similar  to  those  used  In  the  geometrical  optics.  The  paths  of  the  Incident  and 
reflected  waves  are  shown  schematically  In  Figure  42.  Here  0,  H,  PS  and  8  de¬ 
note  the  center  of  the  muzzle  brake.  Its  height,  location  of  the  pressure 
sensor,  and  the  angle  of  Incident  and  reflection  with  respect  to  the  normal 
from  the  ground.  It  must  be  noted  that  the  two  angles  are  identical  which 
guarantees  only  a  single  reflected  wave  can  go  through  the  point  PS.  Following 
the  procedure  commonly  used  In  the  geometrical  optics,  the  reflected  wave  can 
be  regarded  as  another  spherical  wave  generated  by  a  mirror  image  of  the 
muzzle  brake.  Then  the  solution  from  the  method  of  characteristics  Is  used 
simultaneously  for  both  the  Incident  and  reflected  waves  according  to  their 
path  lengths.  Since  the  pressure  Is  a  scalar  quantity,  its  value  at  a  certain 
point  and  at  a  certain  Instant  is  obtained  by  summing  the  contributions  from 
the  two  waves.  Strictly  speaking,  a  superposition  procedure  Is  not  valid  for 
nonlinear  equations,  Equations  (6)  to  (8).  However,  the  difference  of  the 
times  of  arrival  for  the  two  signals  is  small  compared  with  the  time  scale  of 
the  blast  wave.  The  nonlinear  effect  Is  thus  neglected  here  In  obtaining  a 
first  order  solution.  Refinements  will  be  considered  later. 


Figure  42.  Paths  of  Incident  and  Reflected  Waves 
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8.  MATHEMATICAL  MODEL  OF  L'JNG 


8.1  INTRODUCTION 

Gun  blast  Injury  Is  a  complex  phenomenon  and  might  be  very  hazardous  to 
the  crews  who  operate  the  gun.  The  Important  clinical,  pathological  and  physi¬ 
ological  Information,  regarding  the  effects  of  direct  air-blast  Injury  on  the 
biologic  subject,  has  been  extensively  reviewed  (for  Instance,  References  10 
and  11).  In  order  to  quantify  the  damage  to  the  lung  due  to  blast  overpres¬ 
sure,  a  mathematical  model  (Refs.  1  and  12)  was  proposed  to  study  the  body- 
dynamical  and  gas-dynamical  responses  of  the  thoraco-abdoml nal  structures.  In 
the  present  study,  we  utilize  this  model  as  a  mathematical  tool  to  analyze  the 
lung  response  on  the  crew  positions,  resulting  from  the  measured  and  our  cal¬ 
culated  blast  overpressures  respectively. 

8.2  MATHEMATICAL  FORMULATION 

The  lung  Is  simulated  as  a  gas-filled  cavity  whose  volume  represents 
the  total  of  the  alveolar  and  bronchial  gas  (Figure  43).  The  gas  In  the  lung 
Is  assumed  to  benave  as  an  Ideal  gas  and  Is  homogeneous  In  pressure. 
Furthermore,  the  thoracic  and  abdominal  walls  are  assumed  to  be  perfectly 
rigid.  The  effect  on  lung  volume  of  chest-wall  and  diaphragmatic  action  Is 
simulated  by  means  of  two  piston  systems.  Essentially,  the  Inertia,  stiffness, 
tissue  resistance,  and  surface  area  of  the  chest  wall  and  abdomen  are  the  de¬ 
terminate  factors  of  this  action.  In  fact,  these  various  factors  are  Included 
In  the  simulation  of  each  piston  as  a  simple  system  consisting  of  a  mass, 
spring,  damping  friction  and  effective  area.  Apart  from  the  above-mentioned 
assumptions,  linear  response  of  the  spring  Is  assumed  to  render  our  calcula¬ 
tion  within  the  elastic  limit.  The  equations  governing  the  motion  of  the  pis¬ 
tons  are  written  as  follows 

Mc^c  +  Cc*c  +  *c*c  "  Ac(p  *  pr  ’ 

Maxa  +  Caxa  +  Va  *  Aa^p  "  pt^  *  (25) 
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Here  M,  C,  k,  Xi  A  and  p  are  the  mass,  damping  coefficient,  spring  coeffi¬ 
cient,  linear  displacement,  effective  area  and  the  pressure.  The  subscripts  c, 
a  and  t  denote  the  chest-wall,  abdomen  and  lung,  respectively. 


The  gas  flow  through  the  respiratory  passageways  and  openings  Is  ap¬ 
proximated  by  an  orifice  whose  geometric  area  Is  the  effective  area  presented 
by  the  respiratory  openings  to  the  air  flow  (Flgura  43).  The  gas  In  the  lung 
Is  compressed  or  expands  adlabatlcally  in  accordance  with  the  motion  of  the 
pistons  and  It  escapes  from  or  enters  Into  the  lung  cavity  through  the  ori¬ 
fice.  The  rate  of  change  of  lung  pressure  can  be  described  in  the  following 
equation. 


apt  ,  YPi  dV  .  n  * 

cPT  T"3t  +  Sgn(p  "  pt) 


AP  /jP  -  Ptl 


$gn(p  -  PA) 


1fp>P£| 
If  P  <  Pt  ) 


a  »  6  *  10"4 


an  empirical  constant. 


Here  V  denotes  the  gaseous  volume  of  the  lung,  Ar  denotes  the  effective  area 
of  the  orifice,  r  denotes  the  ratio  of  the  specific  heats,  and  V0  denotes  the 
undisturbed  lung  volume.  The  first  term  in  the  right  hand  side  of  Equation 
(26)  represents  the  adiabatic  compression  process  while  the  second  term  Is  ob¬ 
tained  empirically  for  compressible  flow  through  an  orifice.  It  should  be 
noted  that  the  driving  force  of  Equations  (24)  to  (26)  Is  the  difference  be¬ 
tween  the  external  blast  pressure  and  the  lung  pressure  acting  upon  the  cross- 
sectional  area  of  the  orifice  and  the  piston  areas. 
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The  various  coefficients  used  in 

Equations  (24) 

to  (26)  are  obtained 

from  dimensional  analysis  based  on  the 

measured 

values 

from  a  2.2  kg  rabbit 

(Ref.  12).  The 

inherent  assumption  in  the  dimensional  analysis  is  the  similar- 

ity  of  tissue  density,  tissue  elasticity  and  body  shapes  among  all  animal  spe- 

cies.  The  evaluated  results  for  a  180 

lb.  man 

are  summarized  in  the  table 

below. 

Coefficients 

Units 

Value  for  180  lb  Man 

Mc 

lb 

0.8906 

Ma 

lb 

3.7850 

Ac 

ft2 

0.2893  . 

Aa 

ft2 

0.1113 

cc 

Poundal -sec/ft 

2.7820  x  102 

ca 

Poundal -sec/ft 

1.1130  x  102 

Kc 

Poundal /ft 

1.2170  x  103 

Ka 

Poundal /ft 

3.7860  x  103 

vo 

ft3 

2.9180  x  10“2 

Ar 

ft2 

1.9030  x  10‘4 

Subsequently, 

Equations  (24)  to  (26)  are  solved 

numerically  by  a  fourth  order 

Runge-Kutta  method. 
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